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Final Report 

VLBI - UTILIZING THE SPACE SHUTTLE 
NASA Contract NAS-5-25543 
October 12, 1982 

The uses of the Space Shuttle transportation system for orbiting 
Very-Long-Baseline Interferometry (OVLBI) have been examined, both with 
respect to technical feasibility and its scientific possibilities. The 
study consisted of a critical look at the adaptability of current technology 
to an orbiting environment, the suitability of current data reduction 
facilities for the new technique, and a review of the new sciencfe that is 
made possible by using the Space Shuttle as a moving platform for a VLBI 
terminal in sp.ce. The conclusions are positive in all respects: no tech- 
nological deficiencies exist that would need remedy, the data processing 
problem can be handled easily by straightforward adaptations of existing 
systems, and there is a signficant new research frontier to be explored, 
with the Space Shuttle providing the first step. 

The VLBI technique utilizes the great frequency stability of modern 
atomic time standards, the power of integrated circuitry to perform real- 
time signal conditioning, and the ability of magnetic tape recorders to 
provide essentially error-free data recording, all of which combine to permit 
the realization of radio interferometry at arbitrarily large baselines. 

This has permitted ground-based radio observers to achieve angular resolution 
a thousand times greater than that achieved by Optical astronomy, and has 
revealed a wide variety of interesting phenomena. The size of the Earth 
is a fundamental limitation and the use of an orbiting radio telescope allows 
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one to pass beyond this barrier. Chapter II gives an overview of the 
method. 

The specific advantages of VLBI in space are given in Chapter III. 
Simulations of OVLBI with various choices of near-Earth orbits are given 
in Chapter IV. 

A single antenna in near-Earth orbit can be used in conjunction with 
ground-based radio telescopes in the simplest realization of OVLBI con- 
figurations. This is the mode in which the Space Shu<-tle would be used. 

There is, however, a strong impetus to place antennas in much larger orbits, 
and these possibilities are examined in Chapter V. 

A critical examination of OVLBI must examine the effects of orbital 
perturbations, precession, and decay. These problems are examined in Chapter 
VI, in which it is shown that all such perturbations can be corrected for. 

The question of the best orbit for OVLBI applications can be formulated, and 
the results are discussed in Chapter VII. It turns out that the common 
shuttle orbit with an orbital inclination of 57° is close to the optimum 
choice. 

The state of technical readiness of present VLBI data acquisition, 
handling, and processing systems is discussed in some detail in Chapter VIII. 
The conclusions are favorable: the requirements of OVLBI can be met by 

straightforward adaptations of existing systems. 

The state of technical readiness, and the intense scientific interest 
of the potential subjects of research lead to the conclusion that the next 
phase of work should be the development of actual mission concepts. It is 
recommended, therefore, that mission concept studies be carried out for near- 
Earth orbiting VLBI experiments, including shuttle-based, space platform-based, 
and free-flying missions. The relative costs and benefits of the various con- 
cepts can guide the program development. 


II. A BRIEF DESCRIPTION OF VLBI TECHNIQUES 


The technique of very long baseline radio interf ormetry (VLBI) 
affords astronomers their hignest-resolution look, at cosmic phenomena. 

The basic principle of VLBI is the same as that of connected-element 
radio interferometry - the synthesis of a large aperture by the combination 
of the outputs of two or more widely-separated antennas. 

In connected-element interferometry the outputs of the antennas are 
brought together in real time, while in VLBI they are recorded on separate 
magnetic tapes at each station and played back together at a later time. 
This, flexibility leads to vastly increased angular resolution for VLBI. 
Connected-element arrays are limited to tens of kilometers in size by the 
requirements of real-time data and time-standard transmission and by the 
constraints of geography. For example, the longest baseline at the Ver> 
Large Array of- the National Radio Astronomy Observatory (The VLA) is B x 
- 36 km, or : 610,000 wavelengths 3t the typical wavelength of 76 cm. This 
corresponds to resolution of 56 - VB nax - 0.35 arc seconds. On the other 
hand, VLBI antennas may be separated by as much as an earth diameter, al- 
though considerably less is more typical. If we take B - 2R^, - 12700 km 
as an example at 76 cm the baseline is : 220 million wavelenghts, corres- 
ponding to resolution of ' 1 milliarc ;erond , 350 times better than the VLA. 

In each kind of interferometry the basic measurement for a source in 
a direction § is the correlated output of two antennas separated by an 
instantaneous baseline B. This is called the visibility V of the source, 
and is a (complex) function of the source brightness distribution on the 
sky, B(a,5), and of the perpendicular component of the baseline. 
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3^ = 3 - (5.3)3. Ej^is usually resolved into components along the 
(proj ected-on-the-sky) east and north directions; in units of the wave 
length of observation, these components are referred to as u and v. As 
the earth rotates with respect to the source so does the baseline, the 
projected baseline tracing out an elliptical path in (u,v) - space. The 
parameters of the u-v track depend on the length and orientation of B, on 
the coordinates (a, 5) of the source, and on the time- included in the 
observation. 

The relationship between the source brightness distribution B(a,5) 
and the measured visibility V(u,v) is that of a Fourier transformation, 

00 co 

V(u,v) = j dx J dv B(x,y) exp [ 2iTi (ux+vy) ] , 

where x is an eastward and y a northward angular displacement from the 
source center (a, 6). knowledge of the visibiltv at all (u,v) would enable us 
to reconstruct the source brightness from the inverse transform, 

00 CD 

3(x,y) = j du J dv V(u,v) exp [ -2Tfi (ux + vy)]. 

_OJ ”00 

However, the visibility can be measured only at those (u,v) points which 

correspond to possible projected separations of the various antennas 

which are used. In connected-element arrays such as the ’/LA, the N antennas 

are so arranged that the coverage of the (u,v) plant which is afforded by 

N(N - 1) 

the Earth's rotation of the ; j baselines is very dense out to some 

maximum spacing B /A, and is zero bevond. As a result, a near lv-complete 
max * r 


Fourier inversion of finite resolution is possible. The result M(x,v) (the 
"dirty map") is the convolution of the true source distribution B(x,v) with 


Che point response R(x,y) (the "dirty beam") of the observation, which is 
the Fourier inverse of the u-v coverage, i.e., 

M(x,y) * B*R «JJ dx' dy' B(x', y') R(x-x', y-y') 

R(x,y) - /du/dv exp [ — 2iri (uv + vy)] . 

(Sampled points) 


The beam R is a sharply-peaked function of angular width 59 (the resolution) of 


59 = A/B 

max 


» 


which, because of the incompleteness of the u-v coverage, possesses a low, 
non-zero sidelobe level. 

The effect of the beam sidelobes is to mix the true source brightness 
at a given point with a fraction of the source brightness of nearby points. 
Because the Fourier inverse of the beam R has zeros (it is the u-v coverage) , 
a recovery of B using the Fourier convolution theorem is not possible, so an 
approximate deconvolution of the dirty map is performed in the map plane. This 
procedure, called "cleaning" the map, consists of successively subtracting point- 
source responses from the highest points in the dirty map, storing the amplitudes 
of these "components", and then reconstructing the map by convolving these point 
amplitudes with a "clean" beam having the same width as the central peak of the 
dirty beam, but no sidelobes (Hogbom 1974). With a dirty beam of low sidelobe 
level tills non-linear procedure works well, and it is possible to produce a 
"clean map" of a complicated source which has a dynamic range (the ratio of the 
brightest map feature to the weakest believable one) close to the limit imposed 
by the noise in the receivers. The dynamic range typically achieved with the 
VLA, where N = 27 so the u-v plane is very well sampled, is several hundred 



5 


Co several thousand. However, in VLBI as currently practiced, the antennas 
are few in number (N < 7) and somewhat arbitrarily located on the Earth, so 
the (u-v) plane can be sampled only sparsely. The "holes" which result in 
the (u-v) coverage have the consequence of a dirty beam with a high sidelobe 
level. Under these circumstances the map cleaning does not work very well, 
and the dynamic range of the clean map is severely limited, often ten or less. 

As we shall demonstrate in the next section, one of the important advantages 
of a Space-VLBI observation over current Earth-based VLBI is the tremendous 
improvement in the density of the (u-v) coverage which is obtainable for a 
given source, with the result of greatly-reduced sidelobe levels, and therefore 
significant. 1 v improved maps. 

The second difference between connected-element interf eromet-y and VLBI 
concerns the measurement of the phase of a visibility pcinf V(u,v). In 
connected-element interferometry a single clock keeps time for each of the 
antennas, so the phase visibility points for any one baseline is directly 
comparable to that obtained simultaneously on every other baseline. (Of course, 
the actual measured phase depends on time-dependent instrumental constants 
which must be removed by observation of suitable point sources, for which the 
visibility is independent of (u,v), and for which the phase is defined to be 
zero). In contrast, in VLBI each station has a separate clock. These clocks 
are sufficiently stable that when the outputs of two stations are combined, a 
coherent cross-correlation may be obtained for a time long enough to reach a 
suitable signal-to-noise ratio for each separate visibility point V(u,v). However, 
because of longer-term clock drifts, and other instrumental and atmospheric 
effects, it is not possible to directly compare the phases of '’isbility points 
obtained on different baselines at the same time, or on a single baseline 
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at widely separated times. Thus a direct Fourier inversion approach to map 
making is not very useful. The earliest VLB1 observations were therefore con- 
fined to measurement of the amplitude of the visibility, and map making was 
reduced to modeling the amplitude of the visibiltiy by the sum of a small 
number of elliptical gaussians. 

Recently, more sophisticated VLBI data reduction techniques which do 
make use of the measured phases have been developed. When observations on any 
three baselines which form a closed loop (i.e., from the same three stations) 
are combined in a suitable way the net visibility phase around the loop is 
essentially free of instrumental effects. However, the resulting "closure 
phase" does depend on the structure of the source, and provides constraints in 
the source map which are in addition to those due to the amplitude information. 
When a reasonable number of antennas are involved, the closure phases contain 
almost all of the phase information that would be obtainable in connected-element 
interferometry (in an array of N antennas, we get all but 2/N of the phase 
information). VLBI data reduction based on this idea is called "hybrid mapping", 
and has greatly improved the quality of VLBI maps. (Hybrid mapping is closely 
related to the "self-calibration" technique employed in phase-stable inter- 
ferometry, in which an in: tial map of the source is used to generate a model 
visibility, which is in turn compared to the data and thereby used to correct the 
gains of the antennas. This procedure may be iterated several times, and works 
when there are many more baselines than antennas, so that the problem is well- 
constrained.). These techniques should be applicable to space-VLBI, but their 
use in the case of numerous baselines which are not simultaneously present has 
not been explored in detail. Of course, the use of closure phase will require 
the simultaneous use of a minimum of two ground stations in concert with a 


space VLBI terminal. 


III. ADVANTAGES OF A SPACE-VLBI STATION 


We hive studied the advantages of using one or more space VLBI 
stations in a variety of orbits to observe celestial objects at various 
cositions in the sky. In Lhis section we discuss the advantages of such a 
VLBI system compared to an Earth-bound system. One of the limitations of 
ground-based VLBI is the lack of large north-south antenna separations. This 
constraint, which arises because of the predominantly east-west distribution 
of habitable places on the earth, results .n poor north-south resolution for 
celestial objects of low declination. A second limitation is that most of the 
large astronomical antennas are in the northern hemisphere, making the southern 
half of the sky effectively unobserved by VLBI. The use of a single space 
VLBI station immediately removes the first of these restrictions, and also 
opens up to VLBI all of the southern sky that can be seen from a single (or 
perhaps two) ground-based antenna. The existence of at least two large antennas 
in the southern hemisphere (Parkes 64-m and the JPL Deep Space Network 64-m 
in Australia) means that with space VLBI we can study objects anywhere in 
the sky. A third, fundamental, limitation of ground-based VLBI is the small 
number of antenna stations available. This results in a coarse and sporadic 
coverage of the Fourier transform (u-v) space, since the spacing between 
various u-v tracks is dictated by the separations of the various antennas. 

Recent proposals for a dedicated VLBI array of approximately 10 antennas at 
sites chosen to "optimise" the u-v coverage go part of the way towards solving 
this problem, but the results are limited by the small number of stations that 
is economically feasible, and he difficulty of seeing a given source simul- 
taneously from two ground stations which are far apart. 

A single space VLBI station provides a continuously changing baseline 
to each ground station being used. Such a space terminal in a near-earth 
orbit circles the Earth about 8 times in a 12 hour period, and because the 


“round stations are in different positions relative to the source for 
each orbit, a large number of u-v tracks result from each ground station. 

Figure III-l shows the kind of u-v tracks generated by one orbit and by 
several consecutive orbits. As we demonstrate explicitly below, this 
results in u-v coverage which is decidedly superior to that obtainable 
with current ground-based antennas, and comparable in density and superior 
in resolution to that of the proposed dedicated VLBI arrays. 

The fourth fundamental limitation of ground-based VLBI is that the 
longest baseline obtainable is the diameter of the Earth. In practice the 
maximum baseline is somewhat smaller than this because of the actual locations 
of Earth-based stations and of the requirement of simultaneous visibility of 
the source at both ends of a given baseline. The latter problem is exacerbated 
by the location of most large ground-based antennas at latitudes between 30 and 
50 degrees north. The use of a single space station results in substantially 
increased usable baselines even for near-Eaith spacecraft orbits. Larger orbi.ts 
accessible to free-flying spacecraft result in dramatic increases in baseline, 
with concomitant increases in resolution. 

Limits on Resolutions Due to Interstellar Scattering 

The ultimate limit to the resolution of VLBI is set by interstellar 
scattering. Measurements at low frequencies have shown that a point source 
is broadened into a finite disk of size 

n 

9 S - 150 (A/A 0 ) milliarcseconds (mas) , 

where q = a/ (a - 2) , a is a poorly-known parameter, and X Q = 368 cm corresponds 
to a frequency of 81.5 MHz (see Rickett 1977). Studies of the Crab pulsar 
have shown that a > 3, so that n < 3; a nominal value might be n * 2. 
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Figure III-l . Tracks in (u-v) space produced by one and by several 
consecutive orbits of a space VLB station and a single ground station. 



The results below are quite sensitive to the exact value of n (because 
the VLBI frequencies are so far from 81.5 MHz), so we present results for 
n = 2 and n = 3. 

The longest useful baseline is the one for which the angular resolution 
of the experiment , 


is equal to 


9 r = 4 = 2 x 10 3 A /B, mas, 

r B cm km 

the scattering disk 9 . Since 

s 
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s 


1.1 x lCf 3 A 2 mas 
cm 

3.0 x 10 0 A 3 mas 
cm 


n = 2, 


n = 3, 


the maximum useful baseli* . , 


expressed in Earth radii, is about 


B 

max 



300 A -1 

cm 

( n. 

= 2 ) , 

1.1 x 10 5 A' 2 

C n 

= 3 ) • 


cm 


In Table III-l we present the maximum baseline, and resulting maximum 
angular resolution, as a function of observing wavelength, for the cases 
1=2 and n = 3, alorg with the resolution currently obtainable (9 ) from 

CiD 

ground-based VLBI, assuming B = R . For the case n = 3 (a = 3) , it is clear 
that almost unlimited resolution Is theoretically possible - interstellar 
scattering falls off with decreasing wavelength so fast (0 « A 3 ) that at the 

higher frequencies used for VLBI the scattering is negligible for baselines 
which are tens of thousands of Earth radii. However, if n = 2 (a = 4) , the 
scattering disk remains a constraint for baselines of a few hundred Earth 
radii, and resolution would be limited to a few hundred times that currently 
available. It might be added that a determination of 2 can only be made with 
long baselines and short wavelengths, so space VLBI will lead to important new 
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TABLE III-l 



Maximum 

Useful 

Baseline and 

Resulting 

Resolution 


\ (cm) 

50 

18 

6 

2.8 

2.0 

1.35 

0 (mas) 
GB 

17 

5.9 

2.0 

0.92 

0.66 

0.45 

9 (mas) 

2.8 

0.36 

0.040 

0.0086 

0.0044 

0.0020 


6.0 

17. 

50. 

110. 

150. 

220. 

1 0(mas) 

0.38 

0.017 

6.5xl0 _u 

6 . 6xl0~ 5 

2.4xl0~ 5 

7.4 x 10~ 6 

B itc) 

[ max E 

44. 

340. 

3100. 

1.4xl0 4 

2.8xl0 4 

6 .OxlO 4 
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information about the interstellar medium. 


In the next section we begin our detailed examination of the properties 
of space VLBI with a discussion of the spacecraft orbits involved. 

Specification of the Orbit 

The orbit of a satellite about the Earth is specified by six elements, 
three which give the size and shape of the orbit and the satellites' initial 
position in the orbit, and three which specify the orbit's orientation in 
space. The semi-major axis a and eccentricity e specify the size and shape 
of an elliptical orbit. The satellite's position in the orbit at an initial 
time t Q may be specified by the initial mean anomaly 

M = (t - t ), [modulus 2 t], 

o T op 

where T is the orbital period and t is the time of perigee passage. The 
period of the orbit is related to its semi-major axis by Kepler's third law, 

T : 

0 ’ 

where G = 6.67 x 10 5 cgs is Newton's constant of gravitv and M = 5.98 x lO^'g 

E 

is the mass of the Earth. Numerically, 

T(minutes) = 1.659 x iO - ’* a 3/ “ (km) 
a(km) = 331.3 T 2/3 (rain). 

The standard circular orbits available from the Space Shuttle will have 
altitude h in the range 270 < h < 400 km, so since the mean radius of the 
Earth is = 6731 km, these have periods in the range 89.8 < T < 92.4 min. 

An Earth-synchronous orbit of period 24 hours has a = 42200 km, so h = 35800 km. 
A highly eccentric orbit with a=100,000 km and e = 0.90, corresponding to 
altitudes 3630 km < h < 184000 km has T - 5240 minutes = 3^ 15^ 20 . 
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The orientation of the orbit is described by two angles which specify 
the plane in which it lies and by one which specifies a direction in that 
plane. For orbits around the Earth it is convenient to refer the orbital 
plane to the plane containing the Earth's center and equator, and whose 
intersection with the celestial sphere is the celestial equator, and to 
choose the fiducial direction in that plane to be toward the Vernal Equinox 
y (origin of right ascension, a = O^ 1 ) . The longitude of the ascending node (2 
is the angular distance along the celestial equator between y and the ascending 
node of the orbit, N A , the point at which the plane of the orbit intersects the 
equatorial plane and at which the projected motion of the satellite is from 
the southern into the northern hemisphere. The inclination of the orbit i is 
the angle (0 < i < it) between the orbital plane and the equatorial plane; orbits 
are called prograde if i < t and retrograde if i > -y, corresponding to satellite 
revolution in the same direction and the direction opposite the Earth's rotation, 
respectively. 

The sixth orbital element is the longitude of perigee w, which specifies 
the location of the perigee of the orbit with respect to the Celestial Equator. 
It is the angular distance along the projection of the orbit onto the celestial 
sphere from the ascending node to the point of perigee passage, and satisfies 
0 u) <_ 2tt . When u> = 0 or u = it, the major axis of the orbit lies in the 
celestial plane of the celestial equator, while when to = — or m - -y~ the 

minor axis is parallel to this plane. The geometry of orbital orientation 
in space is illustrated in Figure III-2. 

The velocity v of the satellite in its orbit is 



<s 
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Figure III-2 . Geometry of the orbital elements for a spacecraft in 
orbit about the Earth. 
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where r is the length of the radius vector between the instantaneous position 
of the satellite and the center of the Earth. The maximum velocity occurs 
at perigee, where r is a minimum. 


r 


min 


a(l-e) , 


and is 


v 

max 


GM E ) 1 / 2 f l+e ) 1/2 91( ^E) 1/2 ( 1+e j 1/2 km/s, 

a 1-e ' a 1-e 


while the minimum velocity occurs at apogee. 


r 

max 


a ( 1+e ) , 


and is 


(1+e) 

v = v 

min max (1+e) 


7.91 


R E 1/2 1-e 1/2 

(“D ( 1+^ 


km/ s . 


Thus the circular shuttle orbits have v = 7.7 km/s, a large elliptical orbit 

(a = 100,000 km, e = 0.9) has v = 8.7 km/s, v . = 0.5 km/s, and a 

max min 

geostationary orbit has v = 3.1 km/s. By way of comparison, the velocity 
of a ground-based station at latitude $ is that of the rotating earth. 


Vg = 0.46 cos (km/s). 

As we demonstrate below, the fact that the station velocities appropriate 
to space VLBI are one order of magnitude larger than those encountered 
in ground-based observations necessitates only minor changes in the data 


processing, but enormous improvements in the possible u-v coverage. 


IV. SIMULATIONS OF SPAC^-VLBI OBSERVATIONS 


NEAR-EARTH ORBITS 


In order to better understand the properties — advantages and restrictions- 
of VLBI with one or more antennas in orbit we have performed simulations of 
such observations, varying the source position, spacecraft orbital elements, 
and the number and locations of the ground stations involved. In this section 
we present the results for simulations involving a single space-borne antenna 
in near-Earth orbit. 

The dependence of the resulting u-v tracks on the relative source- 
orbit orientation is demonstrated by the examples in Figures IV-la through IV-ld 
(summarized in Table IV- 2). In each case the source is at a * 0* 1 , 5 = +30°; 
referring to Figure III-2, the source has right ascension corresponding to the 
direction of the Vernal Equinox Cy), and is located 30° above the plane of 
the celestial equator. The orbit is chosen to be circular and to have inclin- 
ation i = 60°; the longitude of the ascending node (£2) is varied from 0° to 
270° in steps of 90°. As a result the orbit is "side-on" (Figure IV-la), 
"face-on" (IV-lb), "side-on" (IV-lc), and "back-on" (IV-ld) to the source 
directions, respectively. The u-v tracks which result from a 24-hour obser- 
vation with a single ground antenna demonstrate how these geometrical rela- 
tionships translate into u-v coverage. It is clear that an orbit which is 
exactly "face-on" to a source with coordinates (a, 5) has elements 



Q = a + 90°; 

opt 


such an orbit results in u-v tracks with the maximal extent both N-S and E-W . 


TABLE IV-1 


Four Orbits Projected on the u-v Plane for a Source at a * 0 , 5 * +30 

Space VLB Terminal* 


Figure 

la 

lb 

lc 

Id 


Ground Stations 

i 

2 

Goldstone 

60° 

0° 

Goldstone 

O' 

o 

o 

90° 

Goldstone 

o 

O 

vD 

O 

O 

GO 

r— 4 

Goldstone 

60° 

270° 


*In all cases the space antenna is in a circular orbit at an 
altitude of 370 km. 
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Figure IV-1: Fourier (u-v) coverage and resulting synthesized beam from 

simulations of VLBI observations of a single source (a = 5 » +30°) from 

four different orbits (i = 60°; SI = 0° , 90°, 180°, and 270°). The angular 
scale assumes that the observing wavelength was 18 cm (v = 1.665 GHz). 
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Figure IV-lc 
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On the right hand side of each figure is a "three-dimensional" plot 
of the dirty beam (point response) corresponding to the u-v coverage shown 
on the left. In these plots the x-axis is East-West, the y-axis is North- 
South, and the z-axis is the amplitude of the response. The angular scale for 
the beam plot will depend on the wavelength used in the particular simulations, 
and can easily be scaled for any desired value. For Figures IV- 1 the wave- 
length used was 18 cm; for the rest of the simulations the figure captions 
contain this information. 

A. Simulations for Three Real Sources 

Several simulations of observations of real astronomical sources are 
summarized in Table IV-2. In each case a series of several ground antenna 
combinations has been used, followed by a spacecraft simulation employing 
three ground stations. The spacecraft orbit was fixed , and the results 
illustrate the various morphologies of u-v tracks which can result. 

For the source 3C84, a high-declination radio galaxy, ground-based VLBI 
produces rather good u-v tracks, even with only three US-based stations (Figure 
IV-2a) . Inclusion of a midwest station (North Liberty, Iowa; Figure IV-2b) helps 
to fill in the "holes", and addition of a European station (Bonn; Figure IV-2c) 
approximately doubles the extent of the u-v coverage. However, use of a single 
space-based antenna, illustrated in Figure IV-2d, clearly improves both the extent 
and density of the u-v coverage. 
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TABLE IV- 2 


Three Real Sources Observed with Various Combinations 
of Ground Stations and Space Stations 


Figure 

Name 

Source 

a 


5 


Ground 

Stations* 

Space VLB Terminal** 

i a 

2a 

3C84 

3 h 

16 m 

+41° 

20' 

HSTK, GB, OVRO 

- 

- 

2b 

3C84 

J h 

16 m 

+41° 

20* 

HSTK, GB, OVRO, NLR0 

- 

- 

2c 

3C84 

3 h 

16 m 

+41° 

20' 

HSTK, GB, OVRO, BONN 

- 

- 

2d 

3C84 

3 h 

16 m 

+41° 

20' 

HSTK, GB, OVRO 

57° 

90° 

3a 

3C273 

12 h 

27 m 

+ 

2° 

19' 

HSTK, GB, OVRO 

- 

- 

3b 

3C273 

12 h 

27 m 

+ 

2° 

19' 

HSTK, GB, OVRO, NLRO 

- 

- 

3c 

3C273 

12 h 

27 m 

+ 

2° 

19' 

HSTK, GB, OVRO, BONN 


_ 

3d 

3C273 

h 

12 

B 

rj 

+ 

2° 

19' 

HSTK, GB, OVRO 

57° 

90° 

4a 

W49 

19 h 

B 

oo 

o 

+ 

9° 

01' 

HSTK, GB, OVRO 
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* HSTK = Haystack Observatory, Westford,MA. 

GB = NRAO, Green Bank, W.VA. 

OVRO = Owens Valley Radio Observatory, Big Pine, CA. 

NLRO = North Liberty Radio Observatory, North Liberty, IW. 

BONN = MPIfR, Bonn, West Germany 

**In all cases the space antenna is in a circular orbit at an altitude of 370 km. 



Figure IV-2 . Fourier coverage and resulting synthesized beam from simulated 


VLB observations of the source 3C84. The angular scale assumes that the ob- 
serving wavelength was 18 cm. 
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Figures IV-3 are Che results of simulations for the low-declination 
quasar 3C273. As summarized in Table IV-2, the antennas involved in IV-3a 
through IV-3d are identical to those used for 3C84 above. Examination 

shows the dramatic effect of changing the source declination from 5 - 41° 
to 5 * 2° : ground-ground baselines involving the usual existing stations have 
almost no North-South baseline projection (v) . As a result, the u-v coverage 
degenerates into East-West lines, and there is essentially no North-South (y) 
resolution of the source. This is reflected in the extremely broad beams in 
Figures IV-3a through IV-3c. However, just as in the case with the high- 
declination source 3C84, use of the same single space antenna produces ex- 
cellent u-v coverage (Figure IV-3d) . [In this case, even tnough the source 
right ascension (a * 12* 1 27 m = 187°) corresponds to a nearly "back-on" view 
of the orbit (fi * 90°), the u-v coverage is good because the difference 
between "face-on" and "back-on" disappears for sources at very low declinations.] 
Finally, in Figures IV-4 we present the results of simulations for the 
interstellar maser source W49. This also lies at low declination (6 = 9°), but 
is about 90° away from 3C273 in right ascension. Figures IV-4a through c 
show that, for a ground-based array, poor North-South resolution occurs for 
5*9°, although it is not as poor as for 3C273. Figure IV-4d shows once 
again that, even for low-declination sources, inclusion of a single space- 
borne antenna results in a dramatic improvement in the density and extent of 
the u-v coverage. Comparison of Figures IV-lc (a = 0* 1 and £} = 180°) and 

IV-4d shows that in this case (for W49) we have the same sense of "edge-on" 
orbits. This is in fact the case; W49 has a 3 19^ 08™ 3 287° and ■ 90°, so 
a - n 3 197°, while the example of Figure IV-lc had a - ft 3 -180° 3 180°. 

Further comparison of Figures IV-lc and IV-4d shows how the additional ground 
stations of the latter case result in denser u-v coverage , and thus in a lower 
sidelobe level than in the case of a single ground station. 
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Figure IV-3. Fourier coverage and resulting synthesized beam from simulated 


VLB observations of the source 3C273. The angular scale assumes that the ob- 
serving wavelength was 18 cm. 
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Figure IV-3c 



Figure IV-4 . Fourier coverage and resulting synthesized beam from simulated 


VLB observation -if the source W49. The angular scale assumes that the ob- 
serving wavlenb'.u was 13 cm. 
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Summarizing this section, we conclude that the use of a single space 
antenna, in a high-inclination (i = 57 ) but otherwise arbitrar'.ly oriented 
("fl doesn't matter much") orbit results in dramatic improvement in the North- 
South resolution for low-declination sources, and in some increase the East- 
West resolution, and in a significant improvement in the density of the u-v 
coverage for all declinations. In all cases the simulations ran for one period 
of 24 hours. Had we chosen a longer period, the ground-ground tracks would be 
unchanged. However, because the spacecraft orbital period and a sidereal day 
are surely incommensurate, a second day would produce space-ground u-v tracks 
which differ from those produced the day before. The density of the u-v 
coverage can, in principle, be made as complete as desired in this way. This 

is a significant advantage over ground-based VLBI, and is independent of the 
number of ground antennas involved in a ground-ground synthesis. Orbital pre- 
cession, considered more fully in Ch. VI, is only a small effect in these examples. 

B. Comparison and Combination with a Dedicated Ground Array . 

The idea of an array of about ten telescopes, at carefully chosen locations, 
dedicated to VLBI, U Very Long Baseline Array, or VLBA) , has been discussed by 
groups at NRAO (1981) and at Caltech (Cohen 1980). A proposal to build such a 
VLBA has been made to ths NSF by NRAO. In this section we compare the properties 
of such a dedicated array with those of space VLBI. Because these two approaches 
to improving the current abilities of VLBI turn out to be quite complementary, 
we also discuss the combination of a VLBA with a space-borne antenna. 

In Table IV-3 we summarize the simulations we have performed using a dedi- 
cated VLB array and various space stations. Fictitious sources at three dec- 
linations, corresponding to three equal-area sections of the Northern celestial 
hemisphere (5 = 64° , 30°, and 6°), were assigned a = 12^ — these choices were 
made to facilitate direct comparison withe the Caltech study. Each source 
was "observed" for 24 hours with a set of ten ground stations - - 
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TABLE IV- 3 


Simulations with a Dedicated VLB Ground Array 
and Various Space VLBI Stations 


Figure 

5 h 

Source (a»12 n ) 

Ground Stations* 

Space ' 
i 

VLB Terminal* 

a 

5a 

+64° 

Network #13 

. 

. 

5b 

+64° 

HSTK, GB, OVRO 

57° 

0° 

5c 

+64° 

HSTK, GB, OVRO 

57° 

270° + 

5d 

+64° 

HSTK, GB, OVRO 

26 °t 

270° + 

5e 

+64° 

Network #13 

57° 

0° 

5f 

+64° 

Network #13 

57° 

270° + 

6a 

+30° 

Network #13 

• 


6b 

+30° 

HSTK, GB, OVRO 

57°++ 

0° 

6c 

+30° 

HSTK, GB, OVRO 

57°t+ 

270° + 

6d 

+30° 

Network #13 

57°tt 

0° 

6e 

+30° 

Network #13 

57°tt 

270° + 

7a 

+ 6° 

Network #13 



7b 

+ 6° 

HSTK, GB, OVRO 

57° 

0° 

7c 

+ 6° 

HSTK, GB, OVRO 

57° 

270° + 

7d 

+ 6° 

Network #13 

57° 

0° 

7e 

+ 6° 

Network #13 

57° 

270° + 


*In all cases the space antenna is in a circular orbit at an altitude of 370 km. 

**HSTK = Haystack Observatory. Westford, MA 
G1 = NRAO , Green Bank, W. VA. 

OVRO =• Owens Valley Radio Observatory, Big Pine, CA. 

Network #13 =■ VLB Array # 13 (Cohen 1980). 

+ * optimum value for this parameter, 
tt = near-optimum value. 
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the "Network 13" chosen as roughly optimal in the Caltech study — as well 
as with various combinations of space antenna orbits and ground stations. 

The resulting u-v coverage and synthesized beams are presented in Figures 
IV-5, IV-6, and IV-7. In all cases the observing wavelength was chosen to 
be that of the Caltech study, X 2.8 era. 

The high declination source (a = 12^, 6 * +64°), observed by the dedi- 
cated array "Network 13", results in the u-v coverage and beam shown in 
Figure IV-5a. It is clear thac the array functions exactly as designed — - 
the coverage is very dense, and the North-South and East-West resolutions are 
comparable, as expected for a high-declination source. In Figure IV-5b we 
contrast this with the results of a three-ground-station, single-space-station 
observation, with the orbit chosen to be "edge-on". Comparison with the 
previous Figure shows one advantage of space VLBI - - the maximum baseline is 
about 50% larger than that obtained with a VLBA. In Figure IV-5c and 
IV-5d we change the spacecraft orbit to have, first, the optimal fl, and second, 
the optimal fl and the optimal i. Comparison of these two Figures with the 
previous one shows that a considerable decrease in sidelobe level is possible 
if the spacecraft orbit is chosen "properly". However, comparison with Figure 
IV-5a shows that the density of the u-v coverage afforded by an array of ten 
ground antennas cannot be achieved by three ground antennas and a single space 
antenna (at least not in a single 24 hour observation), 

In order to make use of the large baselines available from space and the 
density of u-v coverage afforded by a dedicated ground array, we have simply 
combined the two. For an "edge-on" orbit, the results are shown in Figure 
IV-5e, and for an optimal J2 (but not optimal i) orbit they are shown in Figure 
IV-5f. The best features of each approach are realized in their combination — the 
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Figure IV-5 . Fourier coverage and synthesized beam from simulated VLB 
observations of a source at a ■ 12^, 6 * +64°. The angular scale corresponds 
to an observing wavelength of 2.8 cm. (.v * 10.665 GHz). 
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Figure IV-5d 






Figure IV-6 . Fourier coverage and synthesized beam from simulated VLB 
observations of a source at a * 12^, <5 = + 30°. The angular scale corresponds 
to an observing wavelength of 2.8 cm. 
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maximum baseline is considerably longer than for the dedicated array alone, 
and the u-v coverage is much denser than for a spacecraft plus a few existing 
ground stations. An observation made with u-v coverage like that in Figure 
IV-5f would be comparable in u-v density to that obtained with the VLA, but 
with the resolution appropriate to an Earth-sized antenna! 

For a source at lower declination (6 = 30°), the dedicated array does 
net function quite as well (Figure IV-6a) . The North-South extent of the 
u-v tracks is lessened (by sin 5) and the length of each u-v track is de- 
creased (by the shortened mutual visibility time of a source of lower 6). 

By comparison, space-borne observations with "edge-on" and "face-on" orbits 
have considerably greater North-South u-v extent and longer tracks, as evi- 
denced in Figures IV-6b and c. Once again, with only a few ground stations 
there are large "holes", and the beam patterns have larger sidelobe levels 
than one would desire. However, combination of the dedicated array and a 
space antenna yields excellent u-v coverage, for either an "edge-on" (Figure 
IV-6d) or "face-on" (Figure IV-6e) orbit. In fact, it works every bit as 
well at 6 = 30° as it does at S = 64°. 

When the source is moved into the third section of the sky, 6=6°, the 
problems with a purely ground-based become severe (Figure IV-7a). Although the 
antenna locations were chosen to yield significant North-South baselines, in 
particular by putting an antenna in Alaska, the v-extent of the coverage is 
quite limited. And at this 5 the lengths of the u-v tracks are highly re- 
stricted by the short mutual visibility times for two antennas witn a large 
East-West separation. Comparison with Figures IV-5a and IV-6a shows that a 
significant degradation in both resoltuion and coverage completeness occurs 
when a VLBA is used at low declination. However, as illustrated in Figures 
IV-7b and IV-7c, the resolution obtained from space VLBI is not highly degraded 
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Figure IV- 7 . Fourier coverage and synthesized beam from simulated VLB 
observations of a source at a = 12^, 6 ■ +6°. The angular scale corresponds 
to an observing wavelength of 2.8 cm. 
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for low 6. This is because the spacecraft orbit is chosen to be highly 
inclined to the celestial equator, yielding a significant North-South 
spacr ntenna — ground-antenna separations even for a source on the celestial 
equator. (In Section VIII below we discuss the choice of a single "best" 
orbit for space VLBI. A tentative answer, i = 57°, is the basis for our 
choice of orbitd inclinf.tion in this section, and is compatible with the 
orbits obtainable with the Shuttle.) Combinations of ground array and space 
antenna, for "edge-on" and "face-on" orbits, result in Figures IV-7d and 7e. 
Especially for the latter orbit, the results of such a combination are nearly 
as good as those fcr sources at high declinations. 

A quantitative comparison of the beams generated by a VLBA alone, by one 
space antenna and three ground stations, and by the combination of a VLBA 
and a space antenna is made in Figures IV-8 through IV-13. The first three 
Figures, IV-8, IV-9, and IV-10 are effectively cross-sections of the beams. 

These are plots of the maximum absolute value of the beam, evaluated in 
ten annuli surrounding the central peak. This effectively displays the 

amplitude of the beam sidelobes as a function of angular distance from the 
center. 

The root-mean-square (i.e., the fluctuation about zero) of the beam as a 
function of radius from the central peak is plotted in Figures IV-11 through 
IV-13. For the region outside the central peak we have used the same annuli 
as above to calculate the fluctuations in the beam tor each of the simulations 
presented in this section. In addition, we nave plotted the half-width at half- 
amplitude of each central peak; both major and minor axes are shown if they 
differ significantly. 

Study of these figures confirms that the combination of a VLBA and a 
space antenna is superior in resolution and in sidelobe level to a VLBA alone, 
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while a space antenna and a few ground stations has superior resolution 


but poorer sidelobe level. The degradation of the performance of a VLBA at 
low source declination may also be seen by comparing the figures. Finally, 
they also demonstrate that an optimum VLBA-plus-space antenna combination does 
not suffer such a fall-off at low source declination. 


Cross-section in annuli of the absolute maxima of the 
beams in Figures IT-5. C'N.tworkW3") , IV-5c ("Shuttle + 3 stations") 
and IV- 5 f ("Shuttle + Network #13"). 
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Figure IV-9. Cross-section in annuli of the absolute maxima of the 


beams in Figure IV-6a ("Network/* 13") , IV-6c ("Shuttle + 3 stations"), 
and IV-6e ("Shuttle + Network //13") . 
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Fi g ure IV- 10. Cross-section in annuli of the absolute 
beams in Figure IV- 7a ("Network //13"), IV-7c ("Shuttle 
and IV- 7 e ("Shuttle + Network 013"). 
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Figure IV-11 . Root-mean-square beam as a function of the angular 
distance from the central peak, for the parts of the beam outside the 
peak for the cases presented in Figures IV-5 (6 =■ 64°). The horizontal 
bars with tick marks show the half-width at half -amplitude of the central 
peaks; where there are two ticks, they correspond to differing major and 

The labels A - F correspond to the Figures IV-5a through IV-5f. 


minor axes. 


RMS 





Figure IV-12 . Root-mean-square beam as a function of the angular distance 
from the central peak, for the parts of the beam outside the peak in the 
cases presented in Figures IV-6 (6 ■ 30°). The horizontal bars with tick 
marks show the half-width at half-amplitude of the central cores; where 
there are two ticks, they correspond to differing major and minor axes. 

The labels A - E correspond to the Figures IV-6a through IV-6e. 
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Figure IV-13 . Root-mean-square beam as a function of the angular 
distance from the central peak, for the parts of the beam outside the 
peak in the cases presented in Figures IV-7 (<5 ■ 6°). The horizontal 
bars with tick marks show the half-width at half-amplitude of the central 
cores; where there are two ticks, they correspond to differing major and 
uinor axes. The labels A - E correspond to the Figure IV-7a through IV-7e. 
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V. SIMULATIONS FOR LARGE ORBITS 


In the previous section we demonstrated the ability of a single space 
antenna in a near-earth orbit to increase the density and extent of the 
u-v coverage of a VLBI observation. The resolution increase of about 50% 
over purely ground-based observations was due to the increased mutual source 
visibility achieved by a space-ground antenna pair. However, in order to 
produce dramati c inct eases in resolution we need to consider spacecraft orbits 
which extend far from the Earth. We saw in 5 III that interstellar scattering 
will limit the longest useful baseline, but that even a pessimistic estimate 
suggests that baselines of several hundred Earth radii will be useful. In 
this section we examine briefly the properties of a deep-space VLB system in- 
volving two or three satellites in highly-elliptic"l orbits. More than one 
satellite is necessary because when the orbit xe is much greater than an 
Earth diameter, all Earth-spacecraft baselines are approximately the same, 
roughly the shape of the orbit itself. Thus filling the u-v plane depends on 
spacecraft-spacecraft baselines. In order co be able to do so for sources in 
any direction, the orbits need to be properly oriented. In what follows we 
study the use of two or three orbits of semi-major axes a = 100,000 km [so the 

longest Earth-spacecraft baseline is B = 2a = 30 R ] , oriented so a the 

max l 

spacecraft-spacecraft baselines have reasonable projections perpendicular to 
any source direction. In all cases the angular scale is appropriate co a 
wavelength of X 18 cm. 

A. Two Large Elliptical Orbits 

Two spacecraft oriented so that their major axes are both in the celestial 
equator and are mutually perpendicular are used in the first simulations. One 
has its minor axis in the plane of the celestial equator, the other has its 
trmur axis perpendicular to this plane. The two semi-major axes are 
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chosen so that the time AT required for 10 orbits of the satellite in the 
smaller orbit corresponds to that for 9 orbits of the other one. Three ground 
stations spaced roughly around the Earth were used, so that there were almost 
always two 3pacecraf t-ground baselines. Sources at seven locations around 
the celestial sphere were "observed" for this time AT. These simulations are 
summarized in Table V-l, and the results presented in Figures V-la through V-lg. 

With only two spacecraft, it is not possible to produce good (u-v) coverage 
for sources at all (a, 5). This is because orbits of high eccentricity, required 
for practical launches from the near-Earth environment, are "long" in only one 
direction. Thus two orbits can span a given plane, but sources in that plane 
will have u or v limited to the minor axis of one of the orbits. This is amply 
illustrated in Figures V-la, Id, and If, where sources in the celestial equator 
are "observed" with two orbits whose major axes lie in that plane. At higher 
declinations, however, the results are much better, and since such orbits are 
the energetically least expensive, a high-resolution VLB based on this kind of 
orbit might well be a good first step away from the Earth. 

B. Simulations with Three Spacecraft 

With the use of three spacecraft, including one whose major axis is highly 
inclined to the celestial equator, it is possible to get extensive (u-v) coverage 
for all source directions (ci, 5). Simulations based or. the three orbits sum- 
marized in Table V-l arc presented in Figures V-2a through 2g . The source posi- 
tions are the same as those used in the two-spacecraft examples. Comparison 
of the several examples shows that a ’’major-axis projection" is available in 
both u and v for all source positions, as expect d. Here we have AT * 8 
largest-orbit periods = 9 middle-orbit periods = 10 smallest-orbit periods, so 
the (u-v) tracks from spacecraft to spacecraft are "closed" and thus repeat 


81 . 


after AT. However, in a real situation, the orbital periods would be 
incommensurate, and observation for extended periods would be able to "fill" 
the (u-v) plane. It should be noted, however, that rhe time scale is set by 
AT = lOO 1 ?!) * 33 days, so such a (u-v) plane filling observations would take 
several months for orbits with a = 100,000 km. 
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TABLE V-l 


Simulations with Two or Three Elliptical Orbits 


Orbit 

a (km) 

e 

i 


U) 

0 

T (davs) 

#1 

100,000 

0.9 

0° 


0° 

0° 

3.64 

// 2 

92,000 

0.9 

yO 

O 

o 


0° 

90° 

3.28 

Figure * 

V-la 

lb lc 

Id 

le 

If 

Is. 


o (hours) 

0 

0 o 

3 

3 

6 

6 


5 (degrees) 

. 0 

45 90 

0 

45 

0 

45 



Orbit 

a (km) 

e 

i 

UJ 


0 

T (davs) 

//I 

100, O0o 

0.9 

0° 

0° 


0° 

3.64 

if 2 

93,200 

0.9 

0° 

0° 


90° 

3.28 

n 

86,200 

0.9 

90 

270° 


0° 

2.92 

Figure * 

V-2a 
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2d 

2 e 

2f 
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Figures V-l . Simulations of VLB observations with two highly-elliptical 
spacecraft orbits. In all cases the orbits are chose described in Table V-l, 
and the seven source directions are spaced around the celestial equator 
and into the Northern celestial hemisphere in such a way as to get a good 
sample of the (.u-v) coverage which is generated. The angular scale assumes 
an observing wavelength of 18 cm (v = 1.665 GHz). 
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Figures V-2 . Simulations of VLB observations with three highly-elliptical 
spacecraft orbits. In all cases the orbits are those described in Table V-l, 
and the seven source directions are spaced around the celestial equator 
and into the Northern celestial hemisphere in such a way as to get a good 
sample of the (u-v) coverage which is generated. The angular scale assumes 
an observing wavelength of 18 cm (v = 1.665 GHz). 
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VI. EFFECTS OF ORBITAL PRECESSION AND DECAY 


In the ideal case of a satellite subject to only the gravitational 
attraction of a spherically-synnnetric planet, the five elements (a,e ,ft,w,i) 
which specify the size, shape, and orientation of its orbit are independent 
of time. However, in the case of a satellite orbiiting the Earth at finite 
altitude, several auditional forces complicate the picture and cause all 
of these elements to change in time. Here we will be concerned with only 
two most important forces, (i) the gravitational forces that arise because 
of the non-sphericity of the Earth and because of the presence of the sun 
and the moon, and (ii) the aerodynamic forces due to the Earth's upper atmo- 
sphere. As we show below, these additional gravitational forces affect ft and u, 
resulting in very useful slow changes in orbital orientation. The atmo- 
spheric drag, on the other hand, results in orbital decay and the eventual 
reentry of near-Earth satellites, and must the-efore be understood and 
mimimized if a long-lasting (> months) mission is anticipated. 

The effect of the non-sphericity of the Earth and the attraction of the sun 
and moon is to produce an effectively non-invers-square-law gravitaional force. 
For near-Earth orbits the effects of the sun and moon are relatively unimportant, 
resulting in regression of the nodes (ft < 0) and precession of the apsides 
(il > 0) at these rates: 


ft 

w 


Moon Sun 

-0® 43/yr 
+0® 54/yr 


-1® 05/vr 
+1° 32/yr 


Although these effects would of course be important for the precise 
spacecraft tracking and location needed for VLBI, they are of little 
importance (or utility) for the overall plan of a near-Earth space VLBI 
mission. The first-order deviation from the sphericity of the Earth, 
described by its quadrapole moment (dimensionless) of 

J 2 = 1.63 x l(f \ 

causes regression of the nodes and precession of the apsides at the rates 

(2 * -4r- J cosi (““) 2 (1-e 2 ) 2 

l 2 a 

and 

u - — — 5 - J, (5cos 2 i - 1) (— ) 2 (1 - e 2 )“ 2 , 

T 2 a 


where R is the Earth's radius and T the satellite's period. Using Kepler's 
third law to eliminate T and inserting numerical values these become 


(2 


- 10.2 (— ) 


7/2 


cosi (l-e -1 ) 


revolutions /year. 


u) = +5.08 ( — ) '^ 2 (5cos 2 i - 1) (1-e 2 ) 2 revolutions/year. 

3 
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For a circular obit uj is irrelevant (and undefined), but the effect of 
regression of the nodes is substantial, because it changes the orientation 
of the orbit with respect to the sky. In Figure VI-1 we have plotted 0 
versus altitude h - a - R for near-Earth circular orbits of various in- 
clinations. Note that orbits with i < 90° regreca, those with i ■ 90° 
have stationary nodes, and those with i > 90 (retrograde orbits) precess. 

The most useful near-Earth orbits for space VLBI have i = 57° (or i ■ 180° 

- 75° = 123°) in order to optimize the resolution of maps made at all possible 
declinations (see § VII). Orbits available directly from the Space Shuttle 
range in inclination from 28° to 57° (for Kennedy launch). The altitudes 
obtainable depend on payload and on the number of OMS units used (one integral 
unit plus zero to three additional OMS units), and range from 'V 280 to ^ 1000 km. 
The maximum payload faJls roughly linearly froir 'v 22000 kg to zero over this 
altitude range, and depends weakly on inclination, lower inclinations enabling 
larger payloads. Over this entire range of orbits the nodal precession ranges 
from -7.8 revolutions/year (h = 270 km, i = 28°) to -3.3 revolutions/year 
(h * 1000 km, i * 57°). Thus any near-Earth mission entails an orbit which 
precesses a full 360° in between 47 and 111 days, and a nominal one does so 
in 81 days. This means that the orbit is oriented more or less face-on to 
every source in the sky several times each year. This is a great advantage 
for space VLBI, for it means that (1) sources anywhere in th6 skv can be 
well observed with any (high-lnc llnation) near-Earth spacecraft orbit, and 
(2) a given source may be observed with the same optimum uv coverage several 
times per year . The second point is of great interest, for some of the most 
fascinating phenomena observed by VLBI involve rapid time variability. 
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Because of the sensitivity of the nodal precession rate to semi-major 
axis, orbits substantially larger than the ones observed above do not offer 
the same advantages. The locus in (i, h) -space of those circular orbits 
which precess (n > 0) or regress (n < 0) times per year, 

- cosl * < io> C1 + r )7/! • 

is plotted in Figure VI-2 for several values of n. The curves for n * t 1 
make a useful demarcation between fast and slow precession. At the in- 
clinations best for VLBI, i = 57° or 123°, such orbits are at an altitude 
of h - 0.63 R = 4030 km. Much larger orbits, such as earth-synchronous 

ones (h * 42200 km) , have much slower precession rates ( = once per 

23 years in the earth-synchronous case). However, there are larger 
orbits which precess at much faster rates due to the effects of lunar 
perturbations. The use of such orbits for space VLBI has been considered 
by the JPL group, and will not be discussed here. 

The effect of the Earth's upper atmosphere on satellite orbits is to 
decrease all three elements a, e, and i, and therefore also the period T. 
Because of the extreme sensitivity of these effects to the upper atmospheric 
density, and of the latter's large and unpredictable variations with slow 
conditions, a priori calculation of satellite lifetimes against changes 
in a, e and i is not likely to be precise. However, the relevant time scale 
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behaves roughly as follows: 


T « 
C 




h /H 
o 



Here T c is the time for decay to re-entry of a satellite in a circular orbit 

of initial altitude h , CM/A) is the mass-to-area ratio of the satellite, 

o 

and the atmosphere is assumed tc have density scale length H and ground-level 
density p Q . If the orbit is initially elliptical, then the first effect will be 
a circularization of the orbit by. decay of the apogee due to energy loss at 
perigee, a process that always takes longer than the subsequent decay of the 
circular orbit that results. In addition, a satellite initially inclined at 
| i | >0 will suffer a decay toward i - 0. 

The lifetime of a satellite depends very sensitively on its initial 
altitude. In the range of Space Shuttle orbits the equivalent atmospheric 
scale height is between 35 and 70 km, so that a change in initial altitude 

i l /I ^ i_ 

of Ah Q results in an increase of at least a factor of e“ in lifetime. 

Thus the difference between h ■ 270 km and h * 400 km is at least a factor 
of 6 in orbital lifetime. The actual lifetimes would depend on the mass-to- 
area ratio of the specific satellite, but it seems likely that an orbit in 
the 400 km range would be suitable for a mission of several months duration. 

A lifetime of years would require either a higher initial orbit or periodic 
visits and orbital adjustments by the Space Shuttle. 
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Figure VI-1 . Race of nodal precession for near-Earth orbits, for 
various orbital inclinations. 
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Figure VI-2. Number of complete nodal precessions per year (n) as a function of 
orbital altitude (h) and inclination (i). 
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h (km) 



Figure Vl-2 





VII. CHOICE OF A "BEST" ORBITAL INCLINATION 


In this section we address the question, "Is there a single optimum 
orbit for a near-Earth orbiting VLBI station?" From the results of §IV 
we know that for a given source position, an orbit that is face-on to a 
given source produces the overall superior synthesized beam. In addition, 
from §VI we know that near-Earth orbits precesses with respect to the stars 
a few full revolutions per year. This means that the optimum ft for a given 
source will occur several times each year, and therefore the specification 
of the initial value ft Q is unimportant for missions lasting a few months or 
more. Thus the orbital inclination i is the only degree of freedom, and the 
question becomes "Is there a single inclination i Q that produces the "best" 
overall results for sources at all possible locations?" As we show below, for 
a very resonable definition of "best overall results" the answer is in the 


affirmative, and the value of i 

opt 


that we obtain is surprisingly convenient. 


For a source at given coordinates (a, 5), a face-on near-Earth circular 


orbit has the elements 


« - 

a 


a + 


IT 

2 . 



The angle between the normal to this best orbit for this source and the normal 
to an orbit at inclination i, at the time when ft * £ , is 
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n = i-l ( 5 = i+ | 5f — j . 


The result of i M, is that the North-South extent of the synthesized beam 
made with such an orbit will increase as the v-orbital projection decreases 


<$0 ns a 1/cos n = l/sin(i + 6). 


On the other hand, as long as SI : the east-west beam width will not 

depend on i. Thus the solid angle of the beam will vary with (i,6) as 


'■ft “ 60 ^,. “ l/sin(i + o) 


If we now define the "quality" of the beam to be inversely proportional to 
its solid angle. 


= l/£fl « sin(i + 6), 


so that is proportional to the number of beam areas in a synthesized map, 
we may find an overall optimum inclination i t by maximizing the 
integrated quality of all the possible maps we could make with respect to 
i. Since the number of sources of declination 6 is proportional to the 
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area of sky at that declination, so that 


dn _ . 

d6 cos 5 ’ 

the "integral quality"of maps made at all possible directions on the sky 
is proportional to 

tt /2 tt / 2 

I j * / dfi Q^(i,5) = J dfi cos 6 sin (i + 5) . 

The condition that there be a single optimum orbital inclination is that 



Using the identity 


0 . 


2 


ULi 

3 12 


< 0 . 


sin(a + b) = sin a cos b + cos a sin b, 


we easily find that 


/V 2 


T/2 


Ij (i) = sin i / cos 2 6 d6 + cos i / sin 5 cos 6 dfi 


TT . , ^ 1 

■£ sin i -f y cos L. 
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Thus the condition of maximum becomes 


3Ij 


3i 


TT 1 

~T~ cos i - -r sin i „ * 0, 

4 opt 2 opt * 


opt 


or 


opt . 


tan' 2 (— j“) 


5 7® 5 2 


Since 


3 1 , 


3 2 i 2 


sin i - 


~2 cos i. 


it is clear that a maximum is achieved, since 


0 < i < j- . 

A most remarkable fact about this result is that such an orbit is just 

the orbit of largest inclination obtainable from a Florida launch of the 

Space Shuttle. In addition, and more importantly, the curve I (i) is quite 

broad about its peak. This is illustrated in Figure VII-1 , where we plot I (i) . 

1 

In fact, any inclination in the range 
32?0< i < 83°5 


■Ii 


3 i 2 


< 0 , 
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has Ij within 10% of its maximum value. Thus the "integrated quality" of 
maps made at all locations on the sky is not very sensitive to orbital 
inclination, as long as it is above about 30° . 

In Figure VII-1 we also plot the ratio of N - S beamwidth for a source 
of declination 5 using an orbit of i to the possible N - S beamwidth 

optj 

(that obtained at i * i^ * ir/2 - 6), 

e ns (i { ,5) si " “ + 6) 


as a function of 5 for several i. In the case l = i „ , = 57?52 the N - S 

optj’ 

IT 2 / 

extent of the beam is never more than (1 + — ) 1 2 - 1.86 times its minimum. 

2 -r 2 V 2 

At the important region near 6 = 0° the ratio is only — (1 + — — • = 1.19. 

Were we to choose i = 45°, the maximum increase in N - S beamwidth is 
1/cos 45° = 7T~ - 1.41 (occuring at 5 = 0 and 5 = 90°). and the integrated 
quality I ( i = 45°) is only 2.4 percent less than the maximum I ( i ■ i t ) . 

We see that excellent results are obtainable from a wide variety of orbital 
inclinations. 

An alternate definition of quality which weights smaller beams more 
heavily, 


Q 2 a (solid angle) ^ « sin^ (i + 6), 
equivalent to the square of the number of picture elements. 
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leads Co 


TT/2 

I(i) ■ / d6 cos6 sin 2 (i + 5); 


changing variables to 


x * i + <5 


and using the identity 


cos(a-b) = cos a cos b + sin a sin b 


we find directly that 


31 

oi 


opt 


i+ TT/2 . 

/ dx <cosi[sinx cos 2 x] + sini[sin 2 x cosx] > 


The condition 9l/5i = 0 can then be reduced to 

opt 



cos i sin 3 i + sin i cos 3 i + cos 4 i - sin 4 i - 0, 
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which is equivalent to 


tan^ (i) - tan^ (i) - tan (i) - 1 - 0 . 


The only physical solution of this quartic equation is 


i = cos 
opt 

2 


-15- V 5 


10 - 


1/2 


= 58® 28, 


which corresponds to a maximum of I 2 . In this case the peak, of I 2 (i) is not 

quite as wide, but I remains within 10% of its peak value for 

38!5 < i < 79°0 (see Figure VII-1). The ratio of beamwidth at inclination 

i to that at the best inclination for a given source is very similar to 
opt 2 

that for the other choice of Q, and has a worst value (at 5 = 90°) of only 

10 1 / 2 10 . i/ 2 

( — ) = 1.90, and an equatorial (5 = 0°) value of ( ) 


5-1/5 5 + /T 

We conclude that any orbit with inclination satisfying 

40° < i < 6Q° , 


with a slight preference for 


i = 57°, 


will yield excellent overall results, and will not severely compromise any 
region of the sky. 


1.18. 
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Figure VII-1 . Integrated "quality" of synthesized maps of sources 
distributed over the sky, as a function of the orbital inclination i 
of the spacecraft VLB terminal. Also shown is the ratio of North-South, 
to East-West beamwidth for synthesis performed for a source of declination 
5 using various satellite inclinations i. 
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Figure VII-1 
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Vin. CONSTRAINTS ON DATA ACQUISITION. HANDLIN G. AND PR nr F SR TNT. 

A. Mark III Svstem 

The Mark III system is a digital data recording and processing system for 
VLBI developed jointly by NASA (Goddard Space Flight Center) , Haystack 
Observatory, NRAO, and MIT for astrophysical, astrometric, and geophysical 
observations. The Mark III system consists of two major components, 1) a set 
of identical data recording stations, one at each observatory in the network, 
and 2) a special purpose processor, located at Haystack Observatory, for 
replaying, correlating, and analyzing data. Since an orbital VLBI station would 
be used in conjunction with the existing ground systems, it is clear that the 
orbital station should be compatible with the standard Mark III system. In 

i 

particular, the output data stream should be Identical in format to that of 
the Mark III stations, in order that it may be replayed by the processor 
without modification or addition of extra hardware. 

• B. Mark III Recording Stations 

The Mark III recording stations each contain an IF distributor, 14 video 
converters, a data formatter, and a 28 track instrumentation tape recorder. At 
each station the IF distributor is suppplied with a stable 5 MHz clock and the 
IF signal from the telescope. The built-in local oscillators in the converters 
have a range of 100-500 MHz in 10 kHz steps. Each converter produces an upper 
and lower sideband output video signal. The video output bandwidth is selec- 
table with 4, 2,1,. 5 .25 and .125 MHz being available. Each video signal is 
clipped and one-bit sampled at the Nyquist frequency, twice the bandwidth. If 
the output voltage at the time of a sample is positive the bit is a "1", if 
negative, a "0". 
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Each of Che 28 digital data streams is then sent to the data formatter. 

The formatter adds parity bits (1 per byte) and sync blocks. The sync blocks 
contain words for identification of the blocks, the time of the first sample of 
the block, and an error checking code for the block. These blocks are 
inserted in each data channel, insuring that each channel can be analyzed 
independently. Each channel of data is recorded on one track of the 28 track 
tape recorder. Because of the addition of a parity bit each byte, the 
recording rate per track is 9/8 the sampling rate. 

C. Mark III Processor 

The data tapes made at each of the recording stations are brought to 
Haystack Observatory for processing. The Mark III processor consists of a 
modular special purpose hardware correlator and a host computer which controls 
the processing. Tapes from three stations can be processed simultaneously. 

The 28 data streams from each of the 3 tapes are routed in such a way that 
each of 84 modules receives a baseline-track pair of data streams. In 
general, the data can be processed at the sane rate it was 
taken. Each module cross-correlates its data streams over relatively short 
(^2s) intervals, known as parameter periods, passing the results for each 
period to a disk file via the host computer. The host computer provides each 
module with the appropriate input parameters for processing each interval, name 
the time delay, delay rate, fringe rate and fringe acceleration. Upon comp- 
letion of the correlation processing, the output records are searched for fring 
by software routines running on the host computer. The correlator processing 
step reduces the data volume by a factor of M.0^. The major requirement for a 


succesful fringe search is That the parameters used by the hardware correlator 
be sufficiently close to the correct values so that the fringes are not lost 
due to less of coherence or to smearing in time. 

Figure VIII-1 Is a block diagram of a correlator module. Each cf cne two 
input data streams Is passed through decoders, which remove sync blocks, time 
tags and parity bits, restoring the data 3tream to radio source data only. 

The time tags are used to properly synchronize the two tape drives. The 
X-data stream then passes through a programmable delay buffer. This buffer, 
using parameters passed from the host computer, makes the final adjustment in 
delay between the two streams to compensate for the geometrical delay, T , 

O 

between the two recording stations. Due to the Earth's rotation, the time 
delay varies during a parameter period. The delay can be incremented or 
decremented by one data sample interval 15 fmes per parameter period, with 
the update intervals specified by the host computer. 

The host computer also supplies the fri;.gt rate, $, and fringe acceleration 
to the module. These are then used to generate 3-level approximations to 
sine and cosine fur.tions with periods equal to the fringe period which each 
multiply the Y data stream. The cosine and sine streams are then cross correlated 
with the delayed X stream at 8 different lags, resulting in 3 complex cross 
correlations. These cross correlations are accumulated for a parameter period 
and the resulting .-urns are output to the host computer for storage in a disk file. 
The fringe phase raodei, $ and $, must be sufficiently accurate that the 
difference between the predicted and actual fringe phase does not change by more 
than a fraction of a rotation during a parameter period. 
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D. Orbital Stations 


An essential requirement for performing useful VLBI observations is 
motion of the receiving antennae relative to one another. In current obser- 
vations, this requirement is achieved by the rotation of the Earth on its 
axis. In VLBI observations involving one or mo*e orbital receiving stations, 
the relative motions will be primarily a result of the orbital motions of the 
spacecraft carrying the stations. In general, this will result in higher 
relative velocities of the antennae, and it must be shown what effect these 
will have on the taking and analyzing of VLBI data. Furthermore, in order 
to analyze the data one must know, a priori , the relative motion of the antennae 
to certain accuracies. It must therefore also be shown that the relative motion 
of the orbital antennae can be determined to within the required tolerances. 
These tolerances, to some degree, will reflect the capabilities of the pro- 
cessing hardware and the amount of time that is considered reasonable to devote 
to analysis of the data. In *u:l-s section we will discuss these issues with 
regard to the current capability of the Mark III VLBI Processor at Haystack. 

The maximum relative velocity of a pair of antennae on the Earth along 
the line of sight to a source depends on antenna location and source direction. 
The vorst case maximum relative velocity is 0.92 km/s for a pair of antennae 
on opposite sides of the Earth, on the equator, observing a source on the 
celestial equator. Relative velocities < 1/3 this value are more typical. The 
maximum relative velocity that can be obtained with one of the antennae in 
orbit will depend on the particular orbit chosen for the spacecraft. Spacecraft 
in circular near near Earth orbits (92 min orbital period) have orbital 
velocities of ^ 7.7 km/s. Spacecraft in synchronous orbits have velocities 
of ^ 3 km/s. (If a highly eccentric orbit is used, the majority of it will be 


r 
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at relatively low velocities. When the spacecraft is near the Larth the 
velocities may be higher, but this portion of the orbit will be less useful 
for VLBI purposes and could be neglected.) The maximum likely relative 
velocity is therefore ^7.7 + .46 = 8 km/s. Again, typical values for the 
projected line of sight velocities will be substantially smaller, and will 
average 'v 15 times those obtained during strictly ground based observations. 

Doppler Shift of Bandpass 

The observing frequency for each channel is selected by the local oscillator 
and frequency synthesizer at each station. If the stations are in relative 
motion along the line of sight then the passband for one station will be Doppler 
shifted relative to the other. It is necessary to keep the Doppler shift 
significantly smaller than the bandwidth of the channel, or the correlation 
between stations will be lost. The Doppler shift is: 

. v 

Av = v 

c o 

where V q is the nominal observing frequency and v is the relative line of 
sight velocity. The Mark III system has six allowable bandwidths, B= 4, 2, 

1, 0.5, 0.25, 0.125 MHz, of which 2 MHz is the most commonly used. The band 
shifts by 


Av = 


V 

c 
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so chat 


Av (in MHz) = 0.8/ACin an) 
max 

Thus at 18 cm there is clearly no problem as ^ v max is only .04 MHz, much 

smaller than the available bandwidths. At 1 cm, however, Av ~ .8 MHz, and 

max 

no longer negligible compared to the largest available bandwidth. This effect 
could be reduced by varying the local oscillators of the orbital station in 
orderto match a nominal RF passband. This would require somewhat more 
complicated hardware in the orbital video converters and increase the analysis 
complexity slightly. It must be emphasized that the sensitivity loss caused 
by this effect can not be recovered in the post-processing stage, although it 
could be calibrated out. 

F • Correlator Parmeters/Restrictions 

The basic parameters which must be supplied to the correlator hardware 
are the geometric delay (tg) , delay rate (ig) , fringe rate ($) and fringe 
acceleration ($). Various hardware limitations in the correlator will place 
restrictions or requirements on each of these quantities. 

De lay . The nominal delay time must be provided to the correlator so the 

two data streams can be aligned within 1/B in ordpr to retain coherence. The 

complex correlation is computed for 8 lags so the nominal delay must be 

accurate within + 4 data samples. One lag corresponds to the data sample 

interval, t s = 1/2B, where B is the bandwidth, so the required delay accuracy 

is + 4/23 = + 1.0 us(2MHz/B) . The uncertainty in delay, Ax is entirely due to the 

uncertainties in the parallel baseline component, A£, ]( , so that At = A£,| /c, and the 
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baseline must be known to an accuracy of 


A/1 |i < 300 (2MHz/B) maters. 


Since the locations of the ground stations ax - * well determined, this requirement 
translates into a requirement on knowledge of the absolute spacecraft location. 
This requirement is independent of observing frequency and the nature of the 
spacecraft orbit (high or low Earth) . 

TM* degree of accuracy should provide no serious difficulty. Spacecraft 
in high orbits will suffer relatively small disturbance forces and their orbits 
should be easy to determine. For spacecraft in low earth orbit, use of STDN 
or TDRSS tracking will be nearly sufficient. The projected errors for STDN 
and rpacelab are at just the required level. Those for TDRSS and Spacelat 
are a factor of 2 too large. At worst this might require an extra pass through 
the data. Once fringes have been found, use of the VLBI system itself with 
bright calibration sources will determine the orbit extremely accurately. 

Delay Rate . As the distance between stations changes the delay must be 
adjusted by the correlator. The adjustment can only be made in units of 
single data intervals, T g . A total of 15 adjustments to the delay + 1 sample 
interval each can be made during a correlator parameter period. This is a 
restriction imposed by the correlator hardware. The delay rate is proportional 
to the relative line of sight velocity. 



The time, T u , for the delay to change by one sample, T , is 


s _ c 
2Bv il 

8 " 
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Since 15 updates plus the initial value are allowed, the update interval 
is 1/16 of a parameter period, Tpp. so that the parameter period must 
satisfy 



8c 

Bv |, 


0.15 ( 


2MHz 

B 


^ 8km/ s j 


sec. • 


The parameter period is normally chosen as large as possible. 2s, 
in the processing of standard VLBI data. The primary effect of using a 
shorter value is to increase the number of parameter periods contained in an 
observation. Slncer there are a fixed number of output words per parameter 
period, this implies that the required disk storage space will increase. The 
worst case is for v^ax, which yields a factor of 13 increased storage. This 
should not be a significant problem since the amount of disk space available 
to the host computer could easily be expanded to more than the ~1 M byte 
currently available. 

A second effect of the shorter parameter period is that the host computer 
must s'" 'ice each correlator module more frequently. The rate at which the 
servicing can be done will depend to some extent on the degree of difficulty 
of the computations used to derive the various model parameters. This could 
require restructuring the system to process less than the full number of t -acks 
simultaneously. Because of possible restrictions in the data rate from an 
orbital station, it is not clear how many channels a space system might have, 
so this may or may not prove to be a concern. 
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Fringe Rate . The fringe rate must be supplied by the host computer to 
each correlator module. In order for the fringe search routines to find 
fringes in longer segments of data, the total phase rotation between the 
model and actual data must be a small fraction £ of a rotation over a parameter 
period. Since the fringe rate is related to the delay rate by 


* * Vo 


r !l 


8 x 10 s 


V| 


Hz, 


a A (cm) y 8 km/s 

the requirement of small phase rotation over a parameter period, 


. T pp * ^ v ll 

5<t>= / p 6$dt = ■ t < e 

o X pp 


means that the difference between true and model parallel velocities must 
satisfy 

Av '| < E - 

'* T pp 


At 18 cm and Tpp = .15, a nominal value of £ » 0.5 requires knowledge of the 
spacecraft velocity to w’ithin ^.6m/s. This requirement can be easily met by 
TDRSS, but misses by about a factor of 2 for STDN. Again, it should not be 
difficult to obtain the orbit to this precision. In particular, for bright 
calibration sources, it will be easy to measure fringe rates to better than 
10 mHz, which corresponds to v ^ 0.2 cm/s it X ■ 18 cm and 0.01 cm/s 
if A ■ 1 cm. This should result in extremely good orbit determinations. 

The maximum fringe rate achievable in the correlator is the recording 
bandwidth divided by K, a software-detectable parameter satisfying K = 16, 8, 
4, 2, 1 only, or 

^correlator ^ 5 3 * »*• 
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Thus for the maximum near-Earth fringe rate K must be chosen to satisfy 

K < 2.5 X (B/2MHz), 
cm 

so for a 2 MHz recording bandwidth K must be chosen according to: 


X > 6 . 4cm 

, _ . > 3 . 2cm 

6.4cm > X 

3.2cm > X > 1.6cm 
1.6cm > X > 0.8cm 
0.8cm > X > 0.4cm 
0.4cm > X > 0 


K * 16 
K = 8 

K * 4 

K = 2 

K =*1 

not possible. 


The fringe-rate resolution of the correlator is 


sample rate 

2 24 K 


0.238 _B 

K 2MHz } 


millihertz , 


so that the largest possible value of K should be chosen. 
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Fringe Acceleration. The fringe rate will not be constant; there is a fringe 


acceleration of 


$ (Hz/s) * — y 



990 (• 


92 min 


2/3 


/\(cm) 


In the absence of a fringe acceleration correction a phase error of 


h $ T 



11.1 T PP 2 

A(cm) ^0.05 s 


rotations 


will occur by the end of one parameter period. This amount is intolerable 
for most VLBI wavelengths. However, the Mark III system has provision for 
fringe acceleration correction made by the correlator hardware. The fringe 
phase $ is updated every K data samples by an amount d<J>, which is a 24 bit 
word that can be between 2 1 and 2 -23 rotations and at either sign. This 
phase update itself can be updated every L data bytes (8L data samples) by 
incrementing or decrementing the least significant bit by one. Here L is 
only integer between zero and one-eighth of the number of data samples in a 
parameter period. The least significant bit in A<J) is 2 ' 3 rotations, so the 
fringe acceleration mode is 


* 


± 2~ 2 3 


rotations 


correlator 


(time for K samples) (time for 8L samples) 


or 


1 $ 


correlator 1 


2.38 x 10 s Hz/s 
K L 


, fs 
' 4 MHz 
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where f s ■ data sample rate » 2 x bandwidth of one track. Combining 
this with the required fringe accelerator for near-Earth space VLBI 
results in a requirement on K and L: 


K L 


24!. X (cm) (^) 2 


Once K is picked to satisfy the size of the fringe rate for a given X 
amd f s , L follows from this equation. For example, at X7.8 cm; where 
K 3 4 is the best choice, we require that the phase update A$ be upated 
every L * 169 bytes, or a total of 443 times in a 0.15s parameter period. 
This is easily accomplished by the current Mark III system. 

Quantization of the phase update A<J> and of its update leads to an 
accumulated phase error over a parameter period. This may be estimated by 
comparing the quantized "stair-step" approximation to the phase update to 
the true required phase changes. The absolute level of the stair-step 
approximation is certainly accurate to within 4 of the update quantization, 
i.e., to within 2 -21 * rotations. The slope of the stair-step is accurate to 
within one step per parameter period. Thus the i c ^ update (at time t) is 
in error by no more than 
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Since this error is added to <}> every K data intervals, or ( — ^ S — ) times 

in a parameter period, the phase error accumulated by the end of the period 

is (f t /K) 
s pp 
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Inserting numerical values yields a phase error or 
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rotations. 


Even for the minimum value of K = 1 this does not present a problem. 
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